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8 ! ABSTRACT 
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We present the results of 10 /im nulling interferometric observations of 13 
Herbig Ae stars using the Magellan I (Baade) and the MMT 6.5 m telescopes. 
A portion of the observations was completed with the adaptive secondary at the 
MMT. We have conclusively spatially resolved 3 of the 13 stars, HD 100546, 
AB Aur, and HD 179218, the latter two recently resolved using adaptive optics 
in combination with nulling interferometry. For the resolved objects we find 
\Q \ that the 10 /mi emitting regions have a spatial extent of 15-30 AU in diameter. 
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We also have some evidence for resolved emission surrounding an additional 2 
stars (V892 Tau and R CrA). For those objects in our study with mid-IR SEDs 
classifications from Meeus et al. (2001), we find that the Group I objects (those 
with constant to increasing mid-IR flux) are more likely to be resolved, within 
our limited sample. This trend is evident in correlations in the inferred disk sizes 
vs. the sub-millimeter SED slope and disk size vs. fractional infrared luminosity 
of the systems. We explore the spatial distribution and orientation of the warm 
dust in the resolved systems and constrain physical models which are consistent 
with their observational signatures. 

Subject headings: instrumentation: adaptive optics, stars: circumstellar matter, 
stars: pre-main sequence, techniques: interferometric 
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Introduction 



Circumstellar disks have now been observed surrounding numerous pre-main-sequence 
(PMS) stars of intermediate mass. Herbig Ae (HAE) stars are of particular interest, as they 
are the evolutionary precursors to stars such as Vega, which are known to harbor solid mate- 
rial in their circumstellar environments. Observations of HAE stars can therefore provide us 
with information regarding planet formation and the evolution of the circumstellar disks as 
the stars evolve onto the main-sequence. High spatial resolution studies at wavelengths near 
10/im are especially useful, as they can probe material inner parts (few AU) of these disks, 
allowing one to directly observe regions where planets could form in the habitable zones of 
such stars. 

Infrared (IR) excess emission from HA E stars was originally e xplained by the presenc e 



of geometrically thin, optically thick disks (ILada Adamj Il992l ; iHillenbrand et al 



a mode l that has been modified to include central star irradia tion (jCalvet et al.lll99ll ). disk 



19921) 



flaring (jChiang fc Goldreichl 119971 ; iKenyon &: Hartmannl 119871 ) , and other structures to ex- 
plain features in the spectral energy distributions (SEP), such as an i nner hole, inner wall 
heating, and self- shadowing (IDullemondl 120021 ; iDullemond et al.ll200ll ). Alternative expla- 



nations for circumstellar emission take the form of dusty envelopes or envelope plus disk 
models ( Miroshnichenko et al.lll999l ; lHartmann. Kenvon, fc Calvetl 19931 ) as well as optically 



thin "haloes" in combination with disks (IVinkovic et al. 



2006 
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A review of obs ervational res u lts pr ior to 2000 can be found in lNatta et al.l (120001 ) . More 
recently, a study by iMeeus et al.l (120011 . hereafter M01) based on ISO spectroscopy of HAE 
stars suggests that they may be classified into two groups based on the shape of the SED, and 
explainable by differences in disk geometry. We compare our observational results to the M01 
classifications in jjH Recent years have also seen the grow th of interferometric observations 
of HAE stars (see lMillan-Gabet et al.l (120011 ) for a review). iLeinert et al.l (12004 ) completed a 



long-baseline, spectrally resolved survey of seven HAE stars using the Very Large Telescope 
Interferometer at a wavelength of 10 /im. Their study found that the 10 fim emission 
regions to be 1-10 AU in size, and determined a correlation between the mid-IR SED slope 
and the physical size of the disk, where objects with larger emitting regions had redder 
SEDs. They have attributed this correlation to differences in disk geometry. Additional 
interferometric observations of HAE stars have been done in the near-IR and include the 



studies of lEisner et al.l (120031 ) and lMillan-Gabet et al.l (120011 ). Eisner et al. used the Palomar 
Testbed Interferometer at 2.2 /im and inferred the spatial structure and orientation of the 
near-IR disks surrounding several stars. In a study using the Infrared Optical Telescope 
Array, Millan-Gabet et al. successfully resolved several near-IR emitting regions surround 
HAE stars, establishing the a size scale for the near-IR region of 0.5-5.9 AU, and finding 
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that the distribution of material surrounding resolved objects (including AB Aur) seem to 
favor a circularly symmetric geometry. We will compare our results to those of previous 
interferometric st udies in £ffl Observat i ons at longer wavelengths include a mill imeter study 



(11995b . 



of HA E stars by Mannings fc Sargentl (119971 ) and observations of AB Aur by iMarsh et al. 



Nulling interferometry is a technique used to study spatially resolved circumstellar ma- 
terial in the presence of unresolved flux, and represents an ideal opportuni t y to o bserve the 



environments of HAE stars. The technique, first proposed by iBracewelll (119781 ) is imple- 
mented by overlapping the pupils of two telescopes (or two subapertures on a single aperture 
telescope) with an optical path length difference of one-half A between the beams. The result 
of such a configuration is a sinusoidal transmission function with the functional form 

T{0) = sin 2 {7rbe/X) (1) 

where b is the baseline of the interferometer, A is the wavelength of observations and 9 is 
the angular distance from the observed object to the pointing center of the interferome- 



ter (IHinz et al.l 120011 ). During observations, the central destructive fringe is placed on the 
unresolved point source. This allows us to detect spatially resolved emission, effectively 
isolating it from the unresolved stellar flux. Thus, nulling interferometry provides the neces- 
sary contrast to observe faint circumstellar material in the presence of a much brighter star. 
In observing pre-main sequence stars at wavelengths near 10 /im, where the circumstellar 
emission dominates the stellar flux (often by more than 2 orders of magnitude), isolating 
the resolved disk emission from the unresolved disk emission gives us valuable information 
about the spatial distribution of circumstellar dust. This technique can detect material as 
close to the star as one-quarter of the fringe spacing where the light is neither suppressed 
nor enhanced. This corresponds to 0.12 arcseconds for the configuration used on the MMT, 
or 12 AU about a star at 100 pc. This is between 2 and 3 times finer spatial resolution than 
the diffraction limit at the wavelengths observed in this study (8 to 13 /mi). 

In more recent observations at the MMT, nulling interferometry was implemented in 
combination with the MMT's adaptive optics (AO) secondary mirror. The addition of AO 
has benefits for nulling. Specifically, wavefront abberrations introduced by the atmosphere 
that might affect the suppression level of the null are significantly reduced or eliminated. 
The stabilization of the incoming wavefront therefore allows us to precisely tune the path 
length between the arms of the interferometer, allowing for the best possible suppression of 
light, hence observations are made more efficient (see §2]). Additionally, a stable wavefront 
and image allows us to integrate for greater periods of time, making fainter resolved flux 
detectable. The unique nature of the deformable secondary also has benefits in mid-infrared 
observations. Since the secondary mirror of the telescope is the deformable mirror, there is 
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no need for an intermediate set of correcting optics between the telescope and science camera. 
This reduces the number of warm surfaces in the optical path, minimizing background and 
maximizing throughput. More technical details regarding the MMT adaptive secondary can 
be found in iBrusa et al.l (120031 ) . 

In §2] we discuss the target sample and observations, as well as the reduction procedure 
for nulling interferometric data. In £J3] we present the results of the observations and discuss 
them in context of models and previous observations in §H 



2. Observations and Data Reduction 



2.1. Target Selection 



Targets for this survey were taken from a catalog of Herbig Ae/Be stars by iThe et al. 



( 11994J ). The objects were chosen to include stars with a spectral type of B8 or later, as HAe 
stars have been found to have a s ignificantly higher incidence of circumstellar disks than 
Herbig Be stars (INatta et al.ll2000l ). All objects (except one, HD 98922) were chosen to be 
within 250 pc, to ensure our ability to spatially resolve a circumstellar disk, should one be 
present. All targets also have a 12 /im flux (IRAS) greater than 10 Jy, to insure sensitivity to 
any resolved emission. A large flux is necessary for non-AO observations, as we are limited 
to short integration times (see §2.4p . The final sample of 13 stars, along with their stellar 
characteristics, is shown in Table [TJ 

Observations of 13 Herbig Ae stars were made in 2001 August and 2002 May at the 6.5 m 
Magellan I (Baade) telescope at the Las Campanas Observatory, Chile and 2002 November, 
2004 January, and 2005 June at the 6.5 m MMT at Mt. Hopkins, Arizona. On these single- 
aperture telescopes, nulling interferometry is implemented by dividing the aperture into two 
identical elliptical subapertures, ea ch 2.5 x 5 m with a baseline of 4 m. The BracewelL 
Infrared Nulling Cryostat (BLINC; iHinz I (120011 )) is the interferometer which re-images the 
primary mirror as two subapertures, and recombines the beams which an appropriate path 
difference. The recomb ined image is fed to the Mid- Infrared Array Camera 3 (MIRAC3; 
Hoffmann et al. I (119981 )). which provides the final stop for the image. Observations were 
taken at 10.3, 11.7, and in a few cases 12.5 /zm, each with a bandwidth of 10%, and a wide 
10.6 /im (N-band) filter with a 50% bandpass. Additionally, observations of HD 100546 taken 
at Magellan I included direct imaging taken with BLINC-MIRAC at longer wavelengths (18.0 
and 24.5 /im). 
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2.2. Non-AO Observations 

Observations taken at Magellan and during the 2002 MMT run were done without 
AO. Without wavefront correction, atmospheric aberrations randomize the path difference 
between the two arms of the interferometer. Thus data must be taken with short frame times 
(50 ms) in order to freeze out the seeing effects. In this case, images are taken contiguously 
in large sets, typically 500 to a set. This number of frames is sufficient to sample the 
phase between the beams, and thus include several frames with a phase difference close 
enought to one-half wave that higher order spatial errors dominate the level of residual flux 
in the focal plane. The images with the best null (destructive interference) and constructive 
interference in each set are selected in order to evaluate the "instrumental null" (see below). 
Ten off-source sky frames, used to subtract out the sky background, are taken after each 
set of frames on-source. For each science object we take several sets of frames, in order to 
effectively evaluate the null. From the lowest and highest fluxes in each set of frames, we 
calculate the instrumental null which is defined as: 

-A FluX nu ll e d/ ' FluX cons tructive (2) 

and expressed as a percent. 

Observations of each science object are followed by 2 to 4 sets of 500 frames of a point 
source (spatially unresolved) calibrator star in order to establish a baseline for goodness 
of null achieved in the most recent science observations. Calibrator nulls generally ranged 
between 5% and 15%. We use the null achieved on the calibrator to calculate the source 
null defined as the difference between the instrumental null on the science object and the 
calibrator: 

S ^science ^calibrator (3) 

This represents the resolved flux as a percentage of the full flux of the target. A source null 
of zero means that the science object is spatially unresolved. 

Additionally, observations are taken with different rotations of the interferometer base- 
line relative to the sky. Typically, two sets of 500 frames are taken at each rotation. This 
allows us to probe for the presence of an elongated or flattened structure (such as an inclined 
disk) if resolved emission is detected, which is accomplished in the following manner. The 
output of the interferometer is a transmission pattern superimposed on the stellar image. 
This transmission pattern has a sinusoidal functional form shown in §1, with the destructive 
interference fringe directly on axis, and aligned perpendicular to the baseline of the interfer- 
ometer. If a flattened extended structure is present surrounding our science target, and the 
major axis of the emission is aligned parallel to the destructive interference fringe, we would 
expect a deeper null (smaller N) than when the fringe and the major axis are orthogonal. 
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The resulting dependence of N vs. the rotation of the interferometer would be sinusoidal 
with a period of 180°, 

S = a + b * sin(PA + 6) (4) 

The parameters a, b, PA, and 9 are determined by the size of the emitting region, its 
inclination, the position angle of its major axis, and rotation of the interferometer baseline, 
respectively. 

2.3. Observations with AO 

Observations taken in 2004 January and 2005 June at the MMT were done with the 
AO secondary. As described in the Introduction, the use of AO significantly increases the 
efficiency of observations. Since, with AO, we are able to precisely tune the path difference 
between the arms of the interferometer, large sets of fast frames to freeze-out atmospheric 
effects are unnecessary. AO also stabilizes the image, generally resulting in smaller errors in 
the null. In this case, we take sets of 10 frames (usually ~ Is integration) of the object tuned 
to destructive interference, followed by 10 frames of the object in constructive interference. 
Observations of the calibrator are taken in the same manner, with off-source sky frames 
taken after each destructive-constructive pair. Calibrator nulls typically ranged between 3% 
and 6%, indicating an improvment in both null depth and stability compared to non-AO 
observations. 

Table [2] summarizes all observations of science targets, taken with and without AO. 

2.4. Nulling Data Reduction 

For non-AO observations, each on-source frame is sky subtracted using a sky frame 
created by median combining the 10 off-source frames taken immediately after the science 
frames. A custom IDL program is used to extract aperture photometry from each on sky 
subtracted frame. The aperture size for each set of frames is chosen by determining the 
radius at which the source flux disappears into the noise. The noise level is assessed in an 
annular region several pixels wide well outside the aperture. For each set of 500 frames, we 
identify the frames with the smallest residual flux (best null) and the brightest constructive 
image. These two frames are used to evaluate the instrumental null. The null for each set 
of calibrator images are evaluated in an identical way. 

For the frames taken with AO, the best destructive frame in each set is determined 
using aperture photometry in a manner identical to the non-AO procedure described above. 
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The instrumental null is determined by using the best nulled image and a median combined 
composite image of the constructive frames. Instrumental nulls for calibrator stars, taken in 
between science objects, were evaluated in an identical manner. 

Source nulls were evaluated for all data by subtracting the instrumental null of the 
calibrator from the instrumental null achieved on the science object. A table of instrumental, 
calibrator, source nulls, and errors is shown in Table [3j The values of the nulls represent an 
average of the data sets taken at each wavelength and position angle, and errors are the 1 a 
dispersion in the values. 

3. Results 

Out of 13 objects observed, we have conclusively resolved three of the objects, HD 
100546 (at 3 different wavelengths and 5 different rotations) and AB Aur (at one wavelength 
and 5 different rotations) and HD 179218 (at one wavelength and 2 rotations). Additionally, 
we have marginally resolved (at about the 2.5 a level) two stars, V892 Tau and R CrA. 
Source nulls for all objects can be found in Table [31 

One may note that source nulls are negative at a significant level (> 2a) for one set 
of data. Negative source nulls are unphysical, hence the appearance here necessitates an 
explanation. The negative value for the set in question (HD 98922, 11.7 jttm) can be at- 
tributed to an inaccurate measurement of the calibrator null taken immediately after the 
science object. Analysis of the calibrator images shows a slightly "dual-peaked" image of 
the star, indicating a slight misalignment in the beams of the interferometer and/or a tem- 
porary degradation in observing conditions (i.e., bad seeing). This results in a poor null 
(abnormally large value) for the calibrator null, hence a negative value for a source null. 
It is important to note that a poor null in the calibrator cannot result in a false positive 
detection of resolved material, and any detections of resolved emission are checked to ensure 
that they are not a result of poor seeing or misalignment. The negative source nulls are 
included here for the sake of completeness in presenting data. Other negative source nulls 
are present in the results, but their significance is marginal, as the values are negative at 
less than a 2 a significance. The large values of calibrator nulls and errors in these cases 
are likely due to short-lived poor seeing affecting a single observation of the calibrator (two 
or more observations of the same calibrator are averaged to produce the quoted value for 
the calibrator null). Again, the inclusion of these results is for the sake of completeness and 
transparency in our methodology. 
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3.1. Simple Disk Models - Assessing Size and Flaring 

In order to infer the size of the emitting region for each resolved object, we use two 
simple disk models: 1) the intensity of the disk is a Gaussian function with the peak at 
the center; and 2) the source of the emission is confined to an annulus of uniform intensity 
around the star. The models are placed at the distance of the object and convolved with 
the transmission pattern of the interferometer, and a theoretical source null is calculated. 
The size of the emitting region in the model is adjusted until the results best match the 
observationally determined source null. Models are fit for each wavelength at which an object 
was observed. For any object which we have data at several rotations of the interferometer 
baseline, we fit a sinusoidal function of the form described in §2.21 to the S vs. PA relation 
determined from our nulling observa tions. The fitting procedure is described in detail in 



Liu et all (120031 ) and lLiu et all (12005( 1 . 



The extent of the 10 fim emitting region in each of the resolved objects is a powerful tool 
for determining the physical structure of warm dust in the system. Several factors contribute 
to the spatial extent of emission. One factor is the degree of disk flaring, thus we explore 
the effect of flaring on the nulling observations. Another factor affecting the extent of the 
emitting region is grain size, as larger blackbody grains would tend to result in an emitting 
region closer to the star, due to the relative efficiency in emission compared to small grains. 
We also examine the effect of grain size on our observations. 

In order to assess the effect that the degree of disk flaring has on the resolvabil- 
ity of a target, we constr uct a simple face-on flared disk model based on the model of 



Chiang Sz Goldreichl (119971 . hereafter CG97). In our model we vary the degree of flaring by 
adjusting the power-law dependence of the flaring term, H/a (Eq. 10, CG97, where H is 
the height of the disk above the midplane and a is the radial separation between the star 
and a point on the midplane of the disk). The exponent of a is varied from (the flat disk 
case, i.e. constant H/a) to 2/7 (the vertical hydrostatic equilibrium case, i.e. H/a ~ a 2 / 7 ). 
Other parameters assumed in the CG97 model include a dust mass of 1% of the gas mass, 
uniformly mixed, a grain size of 0.1 /xm, and a mass density of 2 g cm -3 . The CG97 model 
was used to determine the brightness from a hypothetical disk, using can onical stellar pa- 



rameters (effective temperature, mass, radius and luminosity) (ICoxi 120001 ) appropriate for 
the spectral types of the targets' parent stars, while varying the amount of flaring in the 
disk. This brightness signature is then convolved with the transmission function of the inter- 
ferometer to predict the source nulls expected for varying degrees of flaring. Figure [1] shows 
the predicted values of the source null as a function of the degree of flaring (represented by 
the value of the exponent of a) for a range of spectral types. We see that the model predicts 
very large source nulls (40-75%) for the maximally flared disks in all cases. This holds true 
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also for the outer part of the model (iDullemond et al.ll200ll . Eq.5) which predicts flaring as 



the disk emerges from the shadowed region. The largest source null for our resolved objects 
approaches 40% for HD 100546, a B9 star. Thus, it appears that if one considers flaring as 
the only factor affecting the observed source null, none of the objects display flaring to the 
degree of the hydrostatic equilibrium case. It is a possibility that this is an indication that 
the warm dust may be in a flatter distribution than expected (possibly due to dust settling 
toward the disk midplane). However, one must keep in mind that other factors (such as the 
orientation of the dust relative to the rotation of the interferometer baseline) can have a 
significant effect on the source null. 

Spatial information about the emission regions surrounding these stars is also important 
in breaking degeneracies in interpreting the SEDs with regard to grain size. For example, 
large grains, with sizes similar to the wavelength of emission, at closer separations the the 
star can manifest themselves in the same way as small, ISM sized grains in a more extended 
distribution. In the model based on CG97 described above, the assumed grain size is 0.1 fim, 
or ISM-sized grains. Thus, for resolved objects, if grain sizes are actually larger than the 
assumed size, this would require a greater degree of flaring than the actual values calculated 
with the model, in order to account for the same level of resolved flux. When comparing 
the relative effect that grain size and flaring has on the temperature profile of the dust, we 
see that varying the flaring has a greater effect on the power law dependence of the effective 
temperature vs. separation from the star. Increasing the flaring from a flat distribution to a 
flared distribution (at vertical hydrostatic equilibrium) results in a change from a T ~ r -°- 75 
relation to T ~ r -°- 4 ; an d to a change in source null of 40-70%, depending on the luminosity 
of the star. 

Table H] summarizes the sizes and flaring parameters that best fit the nulling data, with 
the errors representing the error in source null shown in Table [3J For unresolved objects, a 
maximum size for the emitting region is shown, calculated using the error in the source null. 



3.2. Notes on Resolved Objects 



HD 100546 This object is perhaps the oldest star in our sample, with an estimated 
age exceeding 10 Myr (see footnote, Tabled]). Large fractions of cryst alline silicates suggests 



an evolved disk, and evidence of a giant protoplanet in the system (IBouwman et al.l 12003 



Malfait et al. 1998: Waelkens et al 



19961 ). HD 100546 was resolved with non-AO observa- 



tions (both nulling interferometric and direct imaging) at Magellan I. Nulling observations 
at 10.3, 11.7, and 12.5 fim show evidence for an inclined disk with an orientation of 45° 
from face-on and a radius of about 12 AU, with a slightly larger size resulting from the ring 



-10- 



model (ILiu et al.ll2003l ). Direct imaging at 18.0 and 24.5 /im also show resolved emission 
at a greater separation (15-20 AU) from the star, and verify the orientation of the dust 
disk. We find that the relative sizes of the 10 and 20 /im emission to not agree with the 
relation for a continuous flared disk (CG97). Instead, a disk with a large inner 



-0.5 



T 

clearing(r < 10 AU) would result in the 10 /xm emission being detected at a greater sepa- 
ration than expected, explaining the relative sizes of the emission at different wavelengths. 
This inner clearing could be the result of the formation of a giant planet. Further discussion 
of this model can be found in lLiu et al.l (120031 ). Here we also consider the possibility that the 
relative sizes of the emission r egions may be a c o nsequ ence of an inner self-shadowed region 
in the disk. In the models of iDullemond et al.l (120011 ). one sees that for the shadowing to 
extend out to the radii of the detected 10 /im emission (>10 AU), the height of the inner rim 
would have to be enhanced relative to the height expected, due to direct normal incidence 
radiation. If this was the case, the models also predict a significant decrease in the strength 
of the 10/ini emission feature. As HD 100546 shows a strong emission feature, it appears that 
self- shadowing is not a significant effect in determining the relative sizes of the the resolved 
emission regions in this system. The relatively large source nulls of HD 100546 suggestive a 
moderate degree of flaring even with the smallest grain sizes. This result supports the fact 
that self-shadowing is not a factor in the r egions we observe and is also in agreement with 
the SED analysis of iBouwman et al.l (120031 ). which infers that the vertical height of the disk 
must be enhanced at ~ 10 AU. 

AB Aur The age of AB Aur is estimated to be significantly younger than HD 100546, 
with its SED showing no evidence for crystalline silicates (M01). Observations of AB Aur at 
10.3 fim (N-band) with AO were taken at 5 different rotations of the interferometer baseline 
and, like HD 100546, show a variation consistent with the presence of an inclined disk. A 
fit to the data results in a disk radius of 12-15 AU, with the smaller sizes in this range from 
the Gaussian disk model, and the larger sizes a result of fitting t o the ring model. These 
sizes are consistent with a recent mi d-IR study bv lMarinas et al. I J2006jh_ A disk inclination 
of 45° to 65° and a major axis PA of 15° to 45° are also inferred (ILiu et al.l 120051 ). The 
orientation of the disk (inclination, PA) derived from these mid-IR o bservations differ from 



several near-IR and millimeter observations previously completed (jFukagawa et al.l 12004 



Eisner et al.ll2003l : iMillan-Gabet et al.ll200ll ; Mannings & Sargentlll997f ). This suggests that 
the circumstellar environment may be more complex than a simple disk structure. However, 
the sizes of the disk at different wavelengths f rom t his and the aforementioned studies, in 
addition to the mid-IR study of IChen fc Jural (120031 ) . indicate that the wavelength vs. size 
relation is consistent with a te mperature profile of T ~ r~ 05 expected from a flared disk. A 



full discussion can be found in ILiu et al.l (120051 ). From relatively large source nulls and the 
fit to the flaring model, a moderate amount of flaring can be inferred. 
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HD 179218 One of the most distant object in our sample (244 pc; iPerryman et al. 



(119971 )). the SED of this object shows significant levels of crystalline silicates (M01). Initial 



observations without AO failed to resolve this star, resulting in a source null of 3 ± 3% at 



a PA of 162° ( IHinz et al.l 120011 ). Given the greater precision with the use of AO, follow-up 
observations were made in 2005 June. These follow-up observations have positively detected 
resolved emission at a levels from 3 to 7% (0.7 to 1.6 Jy). A plot of source null vs. PA 
is shown in Fig. [2j An average of the source nulls implies a FWHM of 20 ± 4 AU for the 
Gaussian disk model or diameter of 27 ± 5 AU for a ring distribution of dust. The low 
source nulls may also suggest flaring in the dust which is small; significantly less than the 
vertical hydrostatic equilibrium case. The 5* vs. rotation relation does not show significant 
variation, consistent with circular symmetry, though significant inclinations cannot be ruled 
out, as a hypothetical disk with an inclination of 45° results in a variation in the source 
null of about 3%, within the errors of the measurement. The presence of significant silicate 
emission in the ISO spectra of M01 also seems to rule out a large inclination for the object, 
whic h would result in a drop o ff in emission intensity at wavelengths shorter than about 10 



fim (IChiang fc Goldreichlll999l ). However, we do not make any definite conclusions regarding 



the spatial orientation of the dust. 

V892 Tau Observations at 11.7/xm without AO, and at 10.3/mi, using AO, show 
resolved emission at a level of about 3 Jy from this Herbig Ae source. The emission is 
detected at a PA of 164°, but no information can be derived as to the overall orientation of 
the emission, as data was taken at only one rotation. Using an average of the source nulls 
obtained in the two data sets at 11.7 fj,m and assuming a Gaussian intensity distribution for 
the dust, the FWHM is 18-28 AU. The ring distribution yields a diameter of 25-37 AU. At 
10.3 /zm, these sizes are 12-16 and 17-23 AU, respectively. The relative sizes of the emission 
at the two wavelengths is consistent within errors to the expected T ~ r -°- 5 relation for a 
CG97 flared disk, assuming purely thermal emission. Signs of flaring can also be inferred 
from the higher level of resolved emission in the 11.7 /im observations, but, when considering 
flaring as the sole factor affecting the extent of the emission region, suggests that the dust 
lies in a flatter distribution than the vertical hydrostatic equilibrium case, though this may 
also be an indication of larger grains. 

R CrA This object was observed without AO at Magellan I and shows marginal 
evidence for resolved emission. The object was observed at two rotations of the interferometer 
and one of the two rotations yielded a positive detection at the 2a level. The level of resolved 
emission is 8% and suggests a spatial extent for the dust of 15 ±4 AU using a Gaussian dust 
distribution and 20 ± 4 AU using a ring distribution. As with V892 Tau, the source nulls 
are possibly indicative of less flaring than a hydrostatic equilibrium situation would expect, 
or larger grains than ISM sizes. 
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4. Discussion 

4.1. Trends in Resolved Objects 

In selecting our HAE targets, we included objects with a range of spectral types and 
ages, in the hope that any evidence of evolution in the PMS environment, or differences 
in the circumstellar region due to the stellar mass of the parent star could be probed. 
However, we have found that the resolved objects in our sample appear to have a range 
of stellar characteristics. The resolved objects have spectral types ranging from late-B type 
(HD 100546 and HD 179218) to A6 (V892 Tau), and ages ranging from approximately 0.1 
(HD 179218) to > 10 Myr (HD 100546). Therefore, there seems to be no distinguishing 
characteristic between unresolved and resolved sources in terms of age or spectral type. 
A plot of the inferred disk size (see §3.11) vs. the stellar age (for those objects with age 
determinations in the literature) is shown in Figure [3h and shows no obvious trend. This 
could be due to either: 1) errors in age estimates for these stars or 2) the fact that time is 
not the sole or dominant factor affecting disk evolution. 

An analysis of the spectral energy distributions (SED) of the observed stars, however, 
does suggest a trend in the resolved objects. M01 categorizes Herbig stars into two major 
groups: Group I with large amounts of mid-infrared excess, and Group II with moderate 
quantities of mid-IR excess, descending at longer wavelengths. For those objects with clas- 
sifications in the literature, we find that two out of three Group I objects were resolved in 
our initial observations (we do not include HD 179218 as a resolved object here since it was 
unresolved in our initial, non-AO observations), while zero of five Group II objects showed 
resolved emission. This trend is evident in a plot of disk size vs. the sub-millimeter SED 
slope (Fig. [3b), which M01 finds correlates well with the mid-IR SED grouping of the star 
and could be considered a surrogate for the evolutionary state of the disk. Objects with 
steeper sub-mm slopes (index < -3) app ear to have larger , resolvable disks. A calculation 



of the Kendall's r correlation coefficient (IPress et al.lll992l ) for the SED slope vs. disk size 



relation yields a correlation probability of 94%, whereas the stellar age vs. disk size relation 
discussed in the previous paragraph yields a correlation probability of 43%. We also see a 
correlation (87% probability) between the fractional IR luminosity (calculated by M01) and 
the disk size (Fig. [3b). 

M01 attributes the difference in the SED between these two groups to disk geometry, 
with Group I objects displaying a significant amount of flaring outside of the inner disk, 
while disks in a Group II source have less flaring, a result of shieldin g from an optically thic k 



inner disk, or perhaps self-shadowing from a puffed up inner wall (iDullemond et al.ll200ll ). 



It is conceivable that the flaring in the Group I objects results in the disk intercepting more 
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radiation at greater radii, making it easier to spatially resolve the dust disk. M01 correlates 
the amount of mid-IR excess to the scale height of the disk, with greater excesses a result 
of more substantial flaring in the disk, consistent with the observed trend in our survey. In 
contrast, the 10 fim emission in the Group II shadowed disk will be confined to the inner 
few AU, making it more difficul t for resolved e missio n to be detected. Thus our results are 



in good agreement with those of iLeinert et al.l (120041 ) which found that the mid-IR emitting 



regions were larger for the redder, Group I objects. 

Despite this line of reasoning, there is reason to be cautious before attributing the 
characteristics of the SED groupings to an all-encompassing physical model. M01 make the 
assertion that Group II objects show evidence that they are more evolved, due to large grain 
sizes. This would seem to imply that older objects are less likely to be resolved than younger 
ones, if time is the dominant factor in the evolution of these systems. In our sample, we do 
not find any trends in resolved objects with age. This would imply that either the age of the 
system is a poor indicator of the evolutionary state of the system, or that the ages attributed 
to the stars are in error. We can gain insight into this issue by comparing HD 100546 and 
AB Aur, the two most conclusively resolved objects. By comparing observations at different 
wavelengths for each object from this and other studies, we find that the temperature vs. 
radius relation for the stars is dramatically different. The temperature profile suggests that 
AB Aur and HD 100546 do not have similar circumstellar environments despite the similarity 
in the 10 fim resolved emission. The former is consistent with the T ~ r ~ - 5 relation expected 



from a simp le flared di s 



t CG 97 while the latter shows evidence for a inner disk gap (see details 



in §3.21 and iLiu et al.l (120031 )). Thus, the evidence here suggests that our ability to resolve 
an object may be due to different circumstances in each system and cannot necessarily be 
attributed to similar physical models. 



5. Summary and Conclusions 



We have carried out 10 ^m nulling interferometric observations of 13 Herbig Ae stars, 
and reach the following primary conclusions: 

- We have conclusively resolved warm dust surrounding 3 objects, HD 100546, AB Aur, 
and HD 179218, the latter previously unreported. Both HD 100546 and AB Aur show 
significant variation in so urce null vs. PA, wh ich is evidence for an elongated structure 
such as an inclined disk (ILiu et al.l 120031 . 120051 ). HD 179218 was resolved in recent (2005 
June) observations with AO, and preliminary results suggest little variation of null vs. PA, 
consistent with an axisymmetric distribution (such as a face-on disk) for the dust, though a 
significant inclination cannot be ruled out given the errors in the measurement. 
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- We have found evidence for resolved emission around an additional two HAE stars, 
V892 Tau and R CrA. Both the sources show resolved emission at a level of a few percent 
of the unresolved flux at 10 /im. 

- The spatial extent of the emitting region in the resolved systems range from 15-30 AU 
in diameter, assuming two models: a Gaussian disk and a ring. 

- Both SED slope and fractional IR luminosity appear to be good indicators of the spatial 
extent of circumstellar dust. Although our sample size is small, it appears th at M01 Group 



I obje cts are more resolvable than Group II objects, a result consistent with iLeinert et al. 



(120041 ) . This trend is evident in the correlation between disk size vs. sub-mm spectral slope 



as well as disk size vs. fractional IR luminosity. 

- There is a lack of correlation in disk size vs. stellar age, perhaps due to uncertainty 
in age determination, and/or the fact that time is not the sole or dominating factor in disk 
evolution. 

- Using a model based on CG97, we evaluate the effect of disk flaring on the resolvability 
of the objects. We find that a hypothetical object with flaring consistent with vertical 
hydrostatic equilibrium would produce a very large source null in all objects. The fact 
that the source nulls observed in these objects is not that large suggests that flaring does 
not impact the resolvability of the objects as much as expected. In this latter case, an 
alternative explanation for the correlation between the SED groups and resolvability must 
be determined. One possible explanation for this is the grain size. With grain sizes larger 
than ISM sizes, flaring must be enhanced to account for the same amount of resolved flux. 

- Follow up AO observations should be made on all objects not yet observed with AO, 
in order to better constrain the limits, or possibly detect, their spatial extent. Observations 
of Group lb objects, those without strong 10 /zm silicate emission, would also be interest- 
ing to assess any correlation between the presence of these emission lines and the spatial 
resolvability of the objects, and further constrain the distribution and orientation of the 
disks. 
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Table 1. Target List 



IN cLIIlt: 


DpcL. ±ype 


d(pc) 


log Agei^iviyrj 




rtcio. 


HR 5999 


A5-7 


210 ±40 


5.7 ±0.3 




1,3,4,6 


KK Oph 


A6 


160 ± 30 


6.5 ±0.5 


II 


2 


DK Cha 


A 


~ 200 






1,3 


HD 150193 


Al 


150 ± 30 


> 6.3 


II 


1,4,6 


HD 98922 


B9 


> 540 






1,4 


HD 104237 


A4 


116 ±8 


6.3 ±0.1 


II 


1,4,6 


51 Oph 


AO 


131 ± 15 


5.5 ±0.2 


II 


1,2,6 


R CrA 


B8 


~ 130 






3 


AB Aur 


AO 


144 ± 20 


6.4 ±0.2 


I 


3,4,5,6 


V892 Tau 


A6 


150 ±10 


> 7 




3,5 


HD 100546 


B9 


103 ±7 


>7t 


I 


1,2,4,6 


HD 163296 


A3 


122 ± 15 


6.6 ±0.4 


II 


2,6 


HD 179218 


B9 


244 ± 55 


5.0 ±0.6 


I 


2,6 



Note. 



Re feren ces- 1) SIMBAD; 2) iLeinert et all (12004h: 3) 



Hamaguchi et al.1 (120051 ): 4 ) Ivan den Ancker et all (119981 ); 5) iNatta et al. 



2000 



, and re f. ther ein), 6) iPerryman et al.l (119971 ) . 



Meeus et all (120011 ) SED Group 

^HD 100546 a t tribut ed to Lower Centaurus-Crux OB association by 
de Zeeuw et al. Jl999h with an association age of 16 Myr iMamajek et al.1 



(12002). 



Table 2. Summary of Observations 



Date Star A (/mi) # Frames Intgr. /frame (s) PA (or rotation*) (°) Tel. Ref. 



1-2 


HR 5999 


11.7 


1000 


0.5 


-165(rot.) 


Mag. 


1 


3-4 


» 


» 


11 


5) 


110(rot.) 


ii 


1 


5-6 


KK Oph 


11.7 


1000 


0.5 


-150(rot.) 


Mag. 


1 


7-8 


» 


5) 


11 


ii 


120(rot.) 


ii 


1 


9-10 


DK Cha 


11.7 


1000 


0.45 


-170(rot.) 


Mag. 


1 


11-12 


HD 150193 


11.7 


1000 


0.45 


145(rot.) 


Mag. 


1 






10.3 


500 


0.5 


97 


MMT 


2 


13-14 


HD 98922 


11.7 


1000 


0.55 


-40 (rot.) 


Mag. 


1 


15-16 


5) 


10.3 


11 


0.55 


ii 


ii 


1 


17-18 


HD 104237 


11.7 


1000 


0.55 


-60(rot.) 


ii 


1 


19-20 


» 


10.3 


» 


ii 


ii 


ii 


1 


20-21 


51 Oph 


11.7 


1000 


0.11 


135(rot.) 


ii 


1 


22 


)> 




500 


» 


45 (rot.) 


ii 


1 


23-24 


R CrA 


11.7 


1000 


0.11 


0(rot.) 


ii 


1 




» 


» 


11 


11 


30 (rot.) 


ii 


1 


25-27 


AB Aur 


10.3 


60 


1.0 


170 


MMTAO 


3 


28-30 


» 


» 


» 


» 


107 


ii 


3 


31-33 






» 




71 


ii 


3 


34-36 


» 


5) 


5) 


11 


131 


ii 


3 


37-39 


5) 




» 


11 


4 


ii 


3 


40-41 


V892 Tau 


11.7 


1000 


0.5 


164 


MMT 


3 


42-43 




» 




5) 


116 


ii 


3 


44-46 


5) 


10.3 


60 


1.0 


160 


MMTAO 


3 



Table 2 — Continued 



Date Star 


A («m) 


# Frames 


Intgr. /frame (s) 


PA (or rotation*) (°) 


Tel. 


Ref. 


HD 100546 


10.3 


1000 


0.5 


-80(rot.) 


MMT 


4 


» 


ii 


ii 


ii 


-50(rot.) 


11 


4 


5) 


ii 


ii 


ii 


-24(rot.) 


11 


4 


11 


ii 


ii 


ii 


10(rot.) 


11 


4 


11 


ii 


ii 


ii 


40 (rot.) 


11 


4 


11 


11.7 


1000 


0.5 


-80(rot.) 


MMT 


4 


11 


ii 


11 


ii 


-50(rot.) 




4 


11 


ii 


11 


ii 


-24(rot.) 


ii 


4 


11 


ii 


11 


ii 


10(rot.) 


ii 


4 


11 


ii 


11 


ii 


40 (rot.) 


ii 


4 


11 


12.5 


1000 


0.5 


-80(rot.) 


MMT 


4 


11 


ii 


11 


ii 


-50(rot.) 


ii 


4 


11 


ii 


11 


ii 


-24(rot.) 


ii 


4 


11 


ii 


11 


ii 


10(rot.) 


ii 


4 


11 


ii 


11 


ii 


40 (rot.) 


ii 


4 


HD 163296 


10.3 


500 


0.5 


94 


MMT 


2 


ii 


ii 


ii 


ii 


10 


11 


2 


HD 179218 


10.3 


500 


0.5 


162 


MMT 


2 


ii 


ii 


ii 


ii 


87 


ii 


2 


ii 


10.6 


90 


1.0 


50 


MMTAO 


1 


ii 


ii 


ii 


1.0 


60 


MMTAO 


1 


ii 


ii 


60 


1.0 


99 


MMTAO 


1 


ii 


ii 


ii 


1.0 


103 


MMTAO 


1 



Note. — References- 1) this paper; 2) Iffinz et all (120011 ): 3) iLiu et all (12005h : 4) iLiu et all (120031 ) 
*The PA probed by the nulling observation depends on both the rotation of the interferometer baseline 
and the paralactic angle of the object when observed. For objects denoted 'rot.' the paralactic angle was 
not recorded, so the rotation of the interferometer (degrees from an arbitrary position) is noted so the 
reader can estimate the relative value of the PA for observations of each object. 



Table 3. Instrumental and Source Nulls 



Star Instr. Null (%) Cal. Null (%) Source Null (%) PA (or rotation) (°) A (/an) Ref. 



HR 5999 


5.7 ±4.6 


5.6 ±2.5 


0.1 ±5.2 


-165(rot.) 


11.7 


1 


» 


1.1 ± 1.1 


» 


-4.5 ±2.7 


110(rot.) 


ii 


1 


KK Oph 


7.5 ±4.2 


5.6 ±2.5 


1.9 ±4.9 


-150(rot.) 


11.7 


1 


» 


4.3 ±2.9 




-1.3 ±3.8 


120(rot.) 


ii 


1 


DK Cha 


10.3 ±4.2 


8.4 ± 1.5 


1.9 ±4.4 


-170(rot) 


11.7 


1 


HD 150193 


10.6 ±2.6 


6.9 ±0.2 


3.7 ±2.6 


145(rot.) 


11.7 


1 


» 


13 ±5 


13 ±2 


0±5 


97 


10.3 


2 


HD 98922 


5.9 ±3.3 


24.1 ± 7.6 


-18.2 ±8.3 


-40 (rot.) 


11.7 


1 




15.2 ±5.8 


12.3 ±3.9 


2.9 ±7.0 


ii 


10.3 


1 


HD 104237 


9.2 ±6.7 


15.5 ± 1.0 


-6.3 ±6.7 


-60(rot.) 


11.7 


1 


» 


10.3 ± 1.9 


15.4 ±3.2 


-5.1 ±3.7 


ii 


10.3 


1 


51 Oph 


11.9 ±1.1 


22.9 ± 10.4 


-11.0 ± 10.5 


1 


11.7 


1 


)> 


7.3 ±5.0 


ii 


-15.6 ± 11.5 


2 


ii 


1 


R CrA 


19.8 ±5.1 


18.1 ±2.7 


1.7 ±5.8 


1 


11.7 


1 


n 


26.4 ± 1.9 


ii 


8.3 ±3.3 


2 


11.7 


1 


AB Aur 


25.9 ± 1.5 


5.1 ±0.2 


20.8 ± 1.5 


170 


10.3 


3 


» 


18.3 ± 1.4 


ii 


13.2 ± 1.4 


107 


ii 


3 




24.1 ±2.9 


ii 


19.0 ±2.9 


71 


ii 


3 


» 


20.9 ± 1.5 


ii 


15.8 ± 1.5 


131 


ii 


3 




31.7 ±3.1 


ii 


26.6 ±3.1 


4 


ii 


3 


V892 Tau 


24.2 ±2.2 


18.5 ±0.8 


5.7 ±2.3 


164 


10.3 


3 


» 


37.9 ±7.2 


22.0 ± 1.3 


15.9 ±7.3 


116 


11.7 


3 




36.4 ±0.2 


ii 


14.4 ± 1.3 


160 


11.7 


3 



Table 3 — Continued 



Star 


Instr. Null (%) 


Cal. Null (%) 


Source Null (%) 


PA (or rotation) (°) 


A (urn) 


Ref. 


HD 100546 


45.5 ± 1.6 


8.6 ±0.1 


36.9 ±1.6 


-80(rot) 


10.3 


4 


» 


39.9 ±7.0 


ii 


31.3 ±7.0 


-50(rot) 


ii 


4 


5) 


41.9 ±2.8 


ii 


33.3 ±2.8 


-24(rot) 


ii 


4 


11 


29.8 ±8.9 


ii 


21.2 ±8.9 


10(rot) 


ii 


4 


11 


31.6 ±0.6 


12.5 ± 1.8 


19.1 ± 1.9 


40 (rot) 


ii 


4 


11 


40.0 ±2.3 


6.4 ±2.4 


33.6 ±3.3 


-80(rot) 


11.7 


4 


11 


26.1 ±2.7 


ii 


19.7 ±3.6 


-50(rot) 


ii 


4 


11 


37.4 ±2.5 


8.6 ±2.0 


28.8 ±3.2 


-24(rot) 


ii 


4 


11 


31.8 ±2.6 


ii 


23.2 ±3.3 


10(rot) 


ii 


4 


11 


24.3 ±2.0 


9.8 ±2.2 


14.5 ±3.0 


40 (rot) 


ii 


4 


11 


35.5 ±2.8 


5.6 ±1.8 


29.9 ±3.3 


-80(rot) 


12.5 


4 


11 


48.8 ±1.3 


ii 


43.2 ±2.2 


-50(rot) 


ii 


4 


11 


30.9 ±1.7 


6.1 ±0.1 


24.8 ±1.7 


-24(rot) 


ii 


4 


11 


26.4 ±1.6 


ii 


20.3 ± 1.6 


10 (rot) 


ii 


4 


11 


16.3 ±5.9 


10.3 ±2.1 


6.0 ±6.3 


40 (rot) 


ii 


4 


HD 163296 


12 ± 7 


13 ± 2 


-1 ±7 


94 


10.3 


2 


ii 


16 ±2 


ii 


3±3 


10 


ii 


2 


HD 179218 


16 ±2 


13 ±2 


3±3 


162 


10.3 


2 


ii 


14 ±2 


ii 


1 ±3 


87 


ii 


2 


ii 


10.5 ±2.4 


3.7 ±0.6 


6.8 ±2.5 


50 


10.6 


1 


ii 


8.5 ±1.8 


3.6 ±0.8 


4.9 ±2.0 


60 


ii 


1 


ii 


7.2 ± 1.2 


ii 


3.6 ±1.4 


99 


ii 


1 


ii 


11.4 ±2.5 


5.9 ±1.2 


5.5 ±2.8 


103 


ii 


1 



Note. — References- 1) this paper; 2) Iffinz et all (120011 ): 3) iLiu et all (120051 ): 4) iLiu et all (120031 ) 



Table 4. Best Fit Models 



Name 


A (/mi) 


Gaussian FWHM (AU) 


Ring diameter (AU) 


Flaring Exp.* 


Rcf. 


HR 5999 


11.7 


< 19 


< 26 


< 0.06 


1 


KK Oph 


11.7 


< 14 


< 19 


< 0.06 


1 


DK Cha 


11.7 


< 17 


< 23 


< 0.05 


1 


HD 150193 


11.7 


< 12 


< 16 


< 0.01 


2 


HD 98922 


11.7 


< 48 


< 68 


< 0.04 


1 


HD 104237 


11.7 


< 12 


< 16 


< 0.07 


1 


51 Oph 


11.7 


< 17 


< 23 


< 0.08 


1 


R CrA 


11.7 


15 ±4 


20 ±4 


0.02-0.07 


1 


AB Aur 


10.3 


27 ±3 


30 ±3 


0.1-0.15 


3 


V892 Tau 


10.3 


14 ±2 


20 ±3 


0.05-0.15 


3 


n 


11.7 


23 ±5 


31 ±6 


0.05-0.15 


3 


HD 100546 


10.3 


24 ±3 


26 ±3 


0.08-0.18 


4 




11.7 


25 ±3 


27 ±3 


0.07-0.18 


4 




12.5 


30 ±3 


33 ±3 


0.02-0.2 


4 


HD 163296 


11.7 


< 11 


< 15 


< 0.08 


2 


HD 179218 


10.6 


20 ±4 


27 ±5 


0.01-0.05 


1 



Note. — References- 1) this paper; 2) lffinz et all (120011 ); 3) [Liu et all (120051 ): 4) lLiu et all (12003 ) 
* in the range to 2/7; see §3.11 
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Fig. 1. — Predicted source null v s. degree of flaring. The d egree of flaring is represented 
by varying the exponent of a in IChiang fc Goldreichl (119971 . Eq.10) from (the flat disk 
case) to 2/7 (the vertical hydrostatic equilibrium case) Spectral types A6 and B8 are shown; 
intermediate types lie between the two curves. The horizontal dashed line represents the 
median (and mean) source null (20%) for all spatially resolved measurements. 
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HD 179218 Source Null vs. PA 
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Fig. 2. — Source null vs. PA for HD 179218, spatially resolved with AO observations. 
The horizontal dashed line is the mean of the measurements, and expected signature (no 
variation in null vs. PA) from a circularly symmetric source. The dotted line shows the 
expected variation in null for a 20 AU diameter, Gaussian intensity distribution inclined at 
45°. 




Figure 1a: Disk Size vs. Stellar Age 
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Figure 1b: Disk Size vs. Sub-mm SED slope 
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Figure 1c: Disk Size vs. Frac. IR Lum. 




Fig. 3. — (a) Disk size (inferred from fitting a Gaussian intensity distribution to the source 
null) vs. stellar age (for ten objects with age determinations in the literature) for our sample. 
Sizes plotted are derived from the 11.7 /im observations except AB Aur (10.3 jum) and HD 
179218 (10.6 /im). No clear age trend is apparent in the data, (b) Disk size vs. sub-mm 
SED slope (for the six objects with values determined in M01). Objects with steeper SEDs 
in the sub-mm (index < -3.0) correlate with Group I objects and appear to have larger, 



